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0.  Introduction 


The  research  done  under  AFOSR-85-0239  falls  into  three  categories: 

(1)  feedback  stabilization  of  distributed  parameter  control  systems; 

(2)  continuum  dynamics  of  materials  exhibiting  phase  transitions; 

(3)  dynamics  of  reacting  chemical  systems. 

Research  in  the  first  two  areas  has  been  an  ongoing  project  of  the 
investigator  over  the  past  several  years.  Research  in  (3)  is  a  new  endeavor 
by  the  investigator  in  collaboration  with  L.  A.  Segel.  The  results  of  the 
research  have  and  will  appear  in  journal  articles.  In  addition  two  Ph.D. 
theses  have  been  completed  under  the  grant  and  one  is  currently  in  progress. 

In  this  report  a  discussion  will  be  given  as  to  the  nature  of  the 
research  and  results  obtained  in  each  of  the  above  three  categories. 


1.  Feedback  stabilization  of  distributed  Darameter  control  systems. 


In  this  research  the  investigator  concerned  himself  with  finding  feedback 
laws  for  control  systems  governed  by  partial  differential  equations.  In 
particular  those  control  systems  which  give  the  dynamics  of  aero-^elastic 
systems  have  been  of  particular  interest.  To  make  the  issue  more  transparent 
the  following  example  drawn  from  the  work  of  Dr.  J.  Hubbard  at  M.I.T.'s  Draper 
Laboratory  may  be  helpful. 

Consider  the  control  system 


32w  +  34w 
3t2  3x4 


0  for  0  <  x  <  L 


3w 

- —  -  0  for  x 
3x 


3t23x 


♦  f(t) 


,3  ,2 

3  w  3  w 


(1.1) 


Here  v(x,t)  denotes  the  displacement  of  a  beam  and  f(t)  is  an  applied 
scalar  boundary  control,  |f(t)|  <  r,  r  >  0.  In  the  absence  of  the  control 


f  the  bean  oscillates  in  an  almost  periodic  fashion.  The  goal  is  to  find  a 
feedback  control 


f(t)  -  F (w(x,t)) 

where  F  is  some  functional  of  w(x,t)  and  its  derivatives  that  will  drive 
w(x,t)  and  wt(x,t)  to  zero  as  t  ♦  ».  This  control  will  then  stabilize  the 
system  and  eliminate  undesired  oscillations.  We  note  the  control  law  must  be 
nonlinear  due  to  the  imposition  of  the  constraint 

|F(w(x,t) ) |  <  r 

In  his  research  the  investigator  reformulated  this  problem  in  an  abstract 
setting  so  as  to  cover  a  large  class  of  problems  in  a  unified  setting. 
Basically  the  idea  was  to  place  the  problem  in  the  format  of  an  evolution 
equation  in  a  Hilbert  space  so  that  (1.1)  could  be  written  in  the  form 

4^  -  Au  +  Bf  (1.2) 

dt 

At 

where  A  is  the  infinitesimal  generator  of  a  linear  semigroup  e  on  a 
Hilbert  space  H,  f  lies  in  a  second  Hilbert  space  E  and  3  is  a  bounded 
linear  operator  mapping  E  to  H,  u  denotes  the  state  of  the  system.  Our 
constraint  is  now  S f I E  <  r. 

In  this  formulation  the  investigator  was  able  to  use  the  tools  of 
dynamical  system  theory  developed  earlier  by  the  investigator,  J.  Ball,  and 
C.  M.  Dafermos  to  find  a  stabilizing  feedback  control.  The  results  appear  in 
the  paper 

"Feedback  stabilization  of  =  Au  Bf  in  Hilbert  space  when 

dt 

If  I  <  r"  by  M.  Slemrod  submitted  to  Mathematics  of  Control,  Signals, 
and  Systems. 

It  is  also  interesting  to  note  tht  by  working  in  a  general  setting  a 
problem  raised  by  A.  V.  Balakrishan  for  the  feedback  stabilization  of  the  NASA 
Spacecraft  Control  Laboratory  Experiment  (SCOLE)  (see  A.  V.  Balakrishnan,  A. 
Mathematical  Formulation  of  a  Large  Space  Structure  Control  Problem,  24th 
Conference  on  Decision  and  Control,  IEEE,  Dec.  1985)  was  solved  in  this  same 
paper  of  the  investigator. 


Continuum 
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The  research  done  in  the  area  of  continuum  dynamics  of  phase  transitions 
centers  on  understanding  dynamics  of  van  der  Waals  like  materials. 


The  reason  for  interest  in  these  materials  is  that  they  provide  the 
simple  mathematical  models  for  materials  wh  ?h  can  exist  in  two  (or  more) 
phases,  e.g.  for  water:  liquid,  vapor,  solid  phase.  More  exotic  materials 
such  as  plutonium  can  exist  in  even  more  phases  (for  plutonium  seven  phases 
are  possible). 

To  see  analytically  what  the  mathematical  issues  are  recall  that  the  van 
der  Waals  constitutive  equation  for  a  compressible  fluid  is  given  by  (see  e.g 
E.  Fermi,  Thermodynamics,  Dover) 

.  R0  a  .  _  , 

p(w,e).__-_  <2.1 

w 

where  p  is  the  pressure,  6  is  the  absolute  temperature,  w  is  the 
specific  volume,  R,  a,  b  are  positive  constants.  It  is  not  hard  to  show 
that  there  is  a  critical  value  of  0,  the  critical  temperature  0crit»  80 
that  for 

0  > ¥■  <w'9>  <  °  »  <2-2 
cnt  <3* 

and  for  0  < 


(w,0)  <  0, 


b  <  w  <  a. 


B  <  w  < 


where  a  and  8  are  constants  a  <  3,  and 


(w,0)  >  0, 


a  <  w  <  B  , 


(the  spinodal  region) 


The  spinodal  represents  an  anamalous  region  where  the  usual  view  that 
pressure  should  be  monotone  decreasing  in  specific  volume  is  violated. 
Furthermore  substitution  of  this  choice  of  p  into  the  isothermal  balance 
laws  of  mass  and  momentum  (in  Lagrangian  coordinates) 

lr  +  T7P(w)  "  0 

(2 

3w  _  3_  Q 

3t  3x  U 

where  u  is  the  velocity  yields  a  mixed  hyperbolic-elliptic  initial  value 


problem  for  6  <  9crit  i‘e*  (2*3)  ls  hyperbolic  on  b<w<<i,  8<w<*  and 
elliptic  on  a  <  w  <  g. 


The  investigator  has  attempted  to  understand  (2.3)  over  the  past  years 
and  written  numerous  papers  on  the  subject  under  the  present  and  previous 
AFOSR  grants. 

Under  the  present  grant  the  main  work  has  been  on  the  resolution  of  the 
Riemann  initial  value  problem  i.e.  solvability  of  the  one  dimensional  initial 
value  problem  for  a  van  der  Waals  fluid  (2.3)  with  initial  data 

u(x,0)  =  u_  u(x, 0)  -  u+ 

x  <  0  ;  x  >  0  .  (2.4) 

w(x,0)  »  w_  w(x,0)  *  w+ 


Here  u_ ,  w_ ,  u+,  w+  are  constants. 

Motivated  by  previous  work  of  the  investigator  an  attempt  to  solve  (2.3), 
(2.4)  was  made  via  the  method  limiting  artificial  viscosity.  Specifically 
(2.3)  was  regular!  .ed  by  the  "viscous"  system 


+  P, 


etu 


XX 


-  u  ■*  EtW 


XX 


(2.5) 


System  (2.5)  has  the  advantage  that  it  admits  solutions  in  the  similarity 


variable 

5  - 

X 

t' 

i.e.  if  we  set 

U  *  u(  £ 

),  w  *  w(£)  and  substitute 

(2.5)  we 

have 

1 

c 

♦  p’  ■ 

EU'  '  , 

-£w* 

-  u'  - 

ew'  •  , 

and  initial  conditions  (2.4)  become 

u(-«)  ■  u_,  w(-«)  -  w_,  u(+«)  ”  u+,  w(+»)  -  w  +  .  (2.7) 

In  his  paper 

A  limiting  viscosity  approach  to  the  Riemann  problem  for  a 
▼en  der  Waals  fluid,  to  appear  Archive  for  Rational  Mechanics 
and  Analysis 

the  investigator  show  for  initial  data  in  different  phases  e.g.  w_  liquid, 
w+  vapor  (2.6)  always  possesses  a  solution  ue,  we.  Also  except  for  one 
special  case  the  limit  as  e  ♦  0  +  of  (uc,we)  exists  *  (u,w)  and  u,  w 


satisfy  (2.3),  (2.4)  i.e.  the  Riemann  problem  is  solved 


In  line  with  this  area  of  research  the  investigator  has  supervised  two 
Ph.D.  theses  at  R.P.I. 

B.  Cassis,  The  method  of  compensated  compactness  applied  to  a  singular 

perturbed  fourth  order  nonlinear  p.d.e.  and  a  mixed  hyperbolic- 
elliptic  system  of  p.d.e.'s.,  submitted  June,  1986. 


M.  Crinfeld,  Topological  techniques  in  dynamic  phase  transitions,  submitted 
May,  1986. 


A  third  thesis  is  being  written  by  a  student  Luis  Leon. 


3.  Dynamics  of  reacting  chemical  systems. 

This  research  is  a  collaboration  effort  with  Professor  L.  A.  Segel  of  the 
Weizmann  Institute  of  Science,  Rehovot,  Israel.  The  problem  is  to  study  an 
approximation  method,  usually  called  the  quasi-steady  state  assumption  (QSSA) 
or  pseudo-steady-state  hypotheses  (PSSH)  which  arises  in  the  study  of  bio¬ 
chemical  kinetics.  In  its  simplest  form,  the  QSSA  deals  with  kinetics 
described  by  systems  of  ordinary  differential  equations  wherein  after  an 
initial  fast  transient,  one  (or  more)  of  the  dependent  variables  may  be 
regarded  as  in  steady  state  with  respect  to  the  instanteous  values  of  the 
other  dependent  variables. 


■'-* -  ~ I'--1’  V* ? V ^ yy  V  V '..■  X.»v.»  u»  <c> 


The  most  studied  example  of  the  QSSA  concerns  a  biochemical  reaction 
wherein  an  enzyme  (concentration  E)  reacts  reversibly  with  another  chemical 
(concentration  S,  the  substrate)  to  form  an  enzyme-substrate  complex 
(concentration  C).  After  some  manipulations  and  elimination  of  E  the  basic 
problem  becomes  analysis  of  the  nonlinear  system  of  ordinary  differential 
equations 


dS 

dt 


-WC>  +  k-!C 


(3.1) 


dC 

dt 


k1(EQ-C)S 


(k_i  +  k2)C 


S{0)  -  SQ,  C(  0 )  =*  0  f 


(3.2) 


where  k_i,  k^,  k2  Eg  are  positive  constants. 


In  the  QSSA  it  is  argued  that  experimental  measurements  are  generally 

performed  after  a  relatively  short  pre-steady  period  but  before  the  substrate 

value  appreciably  decays.  It  is  in  this  period  that  chemists  argue  that  C 

is  approximately  a  constant  so  that  ~  -  0.  But  if  *  0  then  (3.2)  tells 

dt  dt 

us 


C  -  E„  S/(k  +S) 
0  m 


k 

m 


(3.3) 


and  by  (3.1)  S  is  governed  by  the  estimate 


dS 

dt 


-*2v 

(k  +S ) 

m 


(3.4) 


Approximation  (3.3),  (3.4)  constitute  the  QSSA. 

The  result  of  the  research  on  the  QSSA  was  the  following: 

(1)  L.  A.  Segel  has  given  an  argument  based  on  appropriate  choice  of 
non-dimensional  parameters  to  formally  suggest  when  the  QSSA  holds. 

(2)  M.  Slemrod  has  used  arguments  from  the  qualitative  theory  of 
ordinary  differential  equations  to  prove  Segel ' s  formal  conjecture 
is  rigorously  true. 

The  collaboration  was  a  true  meshing  of  two  different  modes  of  thought  to 
solve  a  classical  problem  in  applied  mathematics.  The  authors  believe  they 
have  the  most  general  result  on  the  validity  of  the  QSSA  to  date.  Furthermore 
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the  authors  result  Is  stated  within  the  utmost  mathematical  rigor. 

A  paper  has  been  prepared  for  submission  to  SIAM  Review  entitled  "A  Study 
of  the  Quasi-steady  State  Assumption",  by  L.  A.  Segel  and  M.  Slemrod. 
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